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Abstract. The long-range spectral density correlations (spectral rigidities A^in) and 
related spectral compressibilities) of the E ^ (bi + b2) Jahn- Teller model are found 
strongly nonuniversal with respect to the Hamiltonian parameters and inhomogeneous 
with respect to the choice of a partial energy segment. However, the partial spectral 
rigidities exhibit common features: an anomalous linear part for small n and a 
saturation for large n. The spectral compressibilities found for the partial spectral 
segments and averaged over a whole relevant part of the spectrum cumulate close to a 
well-defined limit pertaining to the semi-Poisson statistics. This is in accordance with 
similar tendencies revealed in the short-range averaged statistical characteristics of this 
model investigated in our previous paper (Majernikova and Shpyrko Phys. Rev. E73 
057202). These features document an inhomogeneous and nonuniversal weakening 
of level repulsions and nonuniversality of level fluctuations on both long and short 
energy scales. The nonuniversality and inhomogeneity of the statistical characteristics 
correspond to a similar behaviour of the chaoticity parameter (a fraction of the chaotic 
phase space of the trajectories) found for the corresponding semiclassical Hamiltonian. 
We ascribe the nonuniversal and inhomogeneous nonintegrability behaviour to the 
changing degree of the brokening of the rotation symmetry when changing parameters 
of our effectively two-dimensional model. It results in a random distribution of the 
respective localized wave functions at all scales up to the size of an available state 
space. The multifractal behaviour of the wave functions is implied from the analysis 
of their (averaged) fractal dimensions which range up to 1.5±0.1 (for 02). This might 
imply the concept of the chaos-assisted tunnelling between the regions of reduced 
degree of stochasticity through regions of high degree of stochasticity. It supports the 
analogy with the two-dimensional Anderson model with marginal-asymptotically far 
metal-insulator transition. The features found allow us to classify the present model 
as a member of the class with a multifractal eigenfunction statistics characteristic for 
the spectra with weakened level repulsion similar to the Anderson model close to the 
metal-insulator transition. 
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1. Introduction 

Investigations of a wide class of complex Hamiltonians show a coexistence of ordered 
and more disconnected chaotic dynamic regions in their phase space in a semiclassical 
hmit [2] . Respective parameters, fractions of the regular and chaotic regions of the phase 
space were used for the evaluation of statistical characteristics of the related quantum 
spectra. Namely, explicit formulae for the nearest-neighbour level spacing distributions 
[3] and the averaged level number variances or the A3 statistics [1[ |5] were illustrated by 
a model of two coupled quartic oscillators. The level repulsion at small level spacings 
was shown to be reduced by the fractionating of the chaotic region which substantially 
complicated related numerical calculations (level spacing distributions). Discrepancies 
between the numerical data and the analytical formulae for the above mentioned short- 
range and also the middle and long-range fluctuation statistics were ascribed especially 
to quantum effects like the tunnelling between regular regions through chaotic barriers 
and the quantum interference effects. Finally, the chaos-assisted tunneling approach to 
the calculation of the level repulsion implied the correct behaviour of the level spacing 
distributions in the limit of small level separations [6]. This kind of tunnelling was 
shown to be connected with the fractal dimension of a two-dimensional phase space due 
to a hierarchical phase space structure of a chaotic system |7]. 

In this paper, we try to show that one of the systems which exhibit fluctuation 
properties typical for a transition region between order and chaos with varying 
(nonuniversal) degree of nonintegrability is the E ® {bi + 62) Jahn- Teller (JT) model. 
The model is represented by two degenerate electron levels coupled to two phonon 
(vibron) modes via two different interaction strengths. Unitary diagonalization of the 
Hamiltonian in the electron space transforms the system onto two highly nonlinearly 
coupled quantum oscillators in two dimensions. 

In a recent paper [1], we have initiated the study of the statistical evidence of 
quantum chaotic patterns emerging in this model. We have shown that the (short- 
range) statistics of the nearest-neighbour level spacings (NNS) in the range of interaction 
parameters apart from the particular symmetry cases (E(8)e JT and exciton model) is 
nonuniversal and tends to a well-defined limit close to the semi-Poisson law P{S) = 
4S exp{—2S). This intermediate statistics between the Poisson (uncorrelated levels) 
and the Wigner-Dyson distribution of a correlated fully chaotic level system was found 
within the frame of random matrix theory as a critical distribution in the metal-insulator 
transition region in the Anderson model of disorder [HI [9] as well as in several models 
with not fully developed chaotic dynamics by Bogomolny [10]. Namely, the semi-Poisson 
distribution of the level spacings in the plasma model was ascribed to the screening 
(restriction to a finite number of nearest neighbours) of the logarithmic pair interaction 
potential. Statistical methods for energy levels and eigenfunctions of disordered systems 
were reviewed by several authors [TT] . 

In section 2, we shall concentrate on the complementary item of the long-range 
statistics for the excited spectra of the ® (61 + 62) JT model with broken rotation 
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symmetry. It comprises the spectral rigidity measure of correlations of the level density 
on scales large when compared to the mean level spacings. Among possible variants of 
this measure one can cite the A3 statistics of Dyson and Mehta [HI El [121 IE] • Namely, 
the A3 value is defined as a (random) quantity which for a given energy interval [—L, L] 
around the value Eq gives the deviation of the least-square linear fit line from the 
staircase function N{E) of the number of levels with energy below E: 



(here the energy is shifted so that the interval is centered around 0; the minimalization 
is performed with respect to parameters A,B). The average of this measure over 
ensemble A3 is the quantity of the present interest. The A3 measure and the correlation 
function of the level number S2(n) = = (iV^) — (A^)^ (fluctuations of the level 

number in an energy band of a width E, n = (N) is the length of the energy interval 
measured by the mean number of the levels inside) are related straightforward as 
A^^n) = ^ Jq {n^ — 2n'^r + r^) T,2{r)dr [Ej. The use of the Aa-statistics instead of 
the level number fluctuations S2 is justified by the fact that its variance is suppressed 
when compared to that of (S^N) |12j . 

The random matrix theory (RMT) predicts the scaling of this measure as logn 
(or logi^^) for the domain of fully developed chaos, meanwhile for the completely 
uncorrelated sequences of levels a linear scaling ~ n is expected [101 [IS] . 

For the Poisson ensemble of an uncorrelated sequence of levels one has S2(r) = r 
and A3(n) = n/15. The number variance T,2{n) was calculated for a set of ensembles 
within RMT supposed to model the behaviour of Anderson type systems [151 IIS] • Its 
asymptotic form was shown to be similar to the uncorrelated case but with the coefficient 
(level compressibility) x < 1, E2(n) ~ for large n. The value x > at weakened level 
repulsion (compared to the metallic limit) refers to the fractal nature of wave functions. 
Restricting to the nearest neighbour interactions of the logarithmic pair potential in 
the Coulomb plasma model Bogomolny et al [lO] found the value of the compressibility 
X = 1/2 which should hold for the models with the semi- Poisson statistics. 

In section 3, we discuss a completing semi classical concept of the degree of 
chaoticity, i.e., of the phase space fraction occupied by the chaotic trajectories /i for 
different segments of the spectra and interaction parameters. A broad region of /i 
from increasing to maximum values in the low and middle parts of the spectra to the 
decreasing again values in the upper parts of the spectra occurs, strongly depending 
on the parameters and location of the spectral segment. These results confirm the 
nonuniversality and varying nonintegrability properties as obtained in section 2. 

In section 4, an independent quantitative statistical analysis of quantum chaotic 
patterns of long-range type is provided by the fractal properties of the eigenfunctions: 
the set of fractal dimensions Dq is defined by the scaling properties of the inverse 
participation ratios Iq{n) related to the eigensolution of the state n of a Hamiltonian as 




(1) 
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Iq{n) = J2 |^n('")P'^- For the sake of application to the excited states of our electron- 

r 

vibron model (defined in section 2) it is suitable to use the spectral representation where 
the relevant substrate is the space spanned by the vibron (phonon) Fock states n [17] 

/,(n)=^|a„|^ (2) 

i 

Here, Cm = ($i|x„); Xn{Qi,Q2) are the exact (numerically calculated for the present 
model) vibron wavefunctions in the space of coordinates Qi and Q2', the base set $i 
is chosen as the set of vibron Fock states - excited states i of the unperturbed two- 
dimensional harmonic oscillator, namely, i is a compound index denoting the direct 
product of the excited states s and t of the vibrons 1 and 2: \i) = \s)i (S) \t)2 ~ 
bl%t\0)i®\0)2. 

The scaling with the fractal dimension Dg in the spectral representation assumes 
that one explores the probabilities Pi^L^n) = J2seiPs{n) = I]sC*sejn5 where the sum is 
taken over the states inside the cube i of dimension 2 comprising P base Fock states 
and L denotes the number of these cubes {L ~ 1/^^)- Then, 

IqAn)^j:nLin)^L-^^^'^-'\ (3) 

i 

where the information dimension Di is understood as usual in the limit g ^ 1 as the 
scaling factor of exp(— X] -Pj logPj) oc L^^. A common assertion found by Kravtsov 
and Muttalib [18] states that the critical statistics in certain class of systems including 
Anderson model at M-I transition is in an intimate relation with the weakly overlapping 
(implying multifractality) wavefunctions. 

2. A3 statistics as a measure of spectral rigidity of the E®(6i + 62) 
Jahn- Teller model 

The E(8>(6i + 62) JT model is defined by the Hamiltonian 

H = n{b{h + + 1)/ + a{b{ + h)a, - (3{bl + b^)^^ (4) 

described, e.g., in our previous papers [H [T71 [19]. In short, the local spinless double 
degenerate electron level is linearly coupled to two intramolecular vibron (phonon) 
modes of the frequency VL by different coupling constants a ^ (3. The pseudospin 
notation with 2x2 Pauli matrices a^., o"^ and unit matrix I refers to the two- level 
electron system. The operators 6j, b\ satisfy boson commutation rules = 5ij 

and define the vibron coordinates Qi = {b\ + bi), i = 1,2. The interaction term oc a 
removes the degeneracy of the electron levels and the term oc (3 mediates the vibron- 
assisted electron tunnelling between the levels (the class of exciton models differs from 
the present one by the absence of the vibron-2 assistance in the tunnelling term ~ /3). 
Two cases of special symmetry comprise the rotation symmetric E®e JT model (a = /3) 
and the polaron model with /5 = or a = 0. Introducing nonequal coupling constants 
the model above generalizes the common E®e Jahn- Teller model. The interest for such 
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a generalization has various sources, in particular: (a) different coupling constants are 
likely to be caused by a spatial anisotropy in a crystal plane; (b) two vibron modes 
need not have necessary the same frequencies. In this case, the Hamiltonian can be 
appropriately rescaled to the present model with equal frequencies but different coupling 
strengths (in what follows we set f2 = 1); and, last but not least, (c) the generalized JT 
system with broken rotational symmetry presents more rich variety of generic chaotic 
properties than its simpler prototype. As it was shown [TJ [T7] it stands closer to the 
generic models of quantum chaotic behaviour. 

In the following we use the vibron eigenf unctions of the transformed Hamiltonian 

H = uHu-^ [miin] 

H = n{Y^ blh + 1) + a{b\ + 6i) - pf3{bl + b2)Rph (5) 

4=1,2 

exactly diagonalized in the electron subspace by the Fulton-Gouterman (FG) unitary 

operator U = ( ' J , where Rph = exp(z7r6|6i) is the vibron reflection 

operator imposing high nonlinearity in the system and the parity p = ±1. It is easy 
to see that in the FG representation the parities p = ±1 are exactly mirror images of 
each other and can be mapped one onto another by a mere change of the sign in the 
definition of the phonon-2 displacements {Q2 thus the system remains doubly 

degenerated and the parity does not have any impact on the properties of the spectrum; 
in what follows we choose p = +1. 

In the representation of radial coordinates in the plane Qi x Q2, Qi = r cos0, Q2 = 
r sin (j) the transformed Hamiltonian ([5]) yields [17] 

1 d f d\ 1 2 1 
H = r — + -r ■ 

2r dr \ dr J 2 2r'^ dcfP' 

+-\/2Q;r(cos — sin (f)Rph)I + p\/2{a — (3)r sin (j)Rph ■ (6) 

The reflection operator Rph in radial coordinates acts as Rph{r, (j)) f{r, (p) = /(r, vr — </>) 
on some f{r,(f)). An example of the numerical eigenfunction to the transformed 
Hamiltonian ([5]) in the space Qi x Q2 is shown in figure [1] Let us note the apparent 
fractal nature of the space distribution of the state amplitud^. The Hamiltonian ([6]) 
commutes with the operator of the angular momentum J = i^bib^ — &1&2) — crj,/2 if 
a = (3. Thus, the eigenfunctions of the symmetric model can be chosen each to pertain 
to a state with a good quantum number |j| = 1/2, 3/2, . . . [17]. 

The last term in ([6]) breaks the rotational symmetry and involves interactions of 
the basis states with different j. Hence, a wavefunction of the full Hamiltonian is now 
distributed over a range of values of the angular momentum (figure [2]). An analogy 

§ It is to be noted that the energy spectrum of the system is invariant with respect to the interchange 
a /3. But the system itself is not invariant with respect to this interchange. For example, the ground 
state is essentially different in the domains of heavy {a > (3) and light (a < f3) polarons 19J. However, 
in the Fock state representation the components of the corresponding wave vectors differ up to the sign 
change which does not affect the fractal properties of wavefunctions in the spectral representation and, 
hence, the symmetry of exposed results with respect to the said transformation. 
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Figure 1. Example of the wavefunction Xn{Qi, Q2) for a = 3, /? = 2 and n — 200. 

to the two-dimensional Anderson model of disorder can be traced if the base functions 
with definite j are considered as 'pseudosites' over which the wavefunction is spread. 
Similarly as in the Anderson model, the properties of the model are determined by 
the relative strengths of the intersite (different j) and the onsite interactions. The 
microscopic reason for this analogy is the assistance of phonons-2 in transitions between 
the levels ([6]) which allows for the emerging of the pseudolattice. To our opinion, it 
is just these transitions with changing symmetry {Rph) which ensured the similarity of 
the results for NNS distribution [T] to the Anderson model close to the M-I transition 
(marginal in two-dimensions). Traces of quantum chaos cause the sequence of these 
coefficients to be distributed in a random fashion with varying n thus supporting the 
analogy to the Anderson model with a priori given random coefficients. 

Calculating the angular part of the matrix elements of the symmetry breaking 
term in ([6]) between the states of harmonic oscillator the selection rules can be obtained 
which assert that the element with definite j is directly connected to the elements with 
j' = j + 2, j — 2 (cf.especially figure 2 with n = 312 where only every odd 'pseudosite' j 
is populated). This reminds us of an analogy with the 'daisy models' of random matrix 
ensembles [20] with dropping every second term which were shown to lead directly to the 
semi-Poisson distribution of NNS. The natural 'length' of our pseudolattice (for a given 
interval of energy) is thus the maximal number of allowed values of angular momentum 
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Figure 2. Examples of distributions of wavefunctions over j-space for a — 2,j3 
Evidently, the states exhibit different extent of localization. 



for a given n (level number), that is L ~ ^/n (a state with main quantum number 
rir = 0, 1, . . . of a two-dimensional harmonic oscillator is ra^ + 1 times degenerated with 
auxiliary quantum number ranging between and n^, hence n ~ nf.). Thus, the analogy 
to the Anderson type models with M-I transition whose intrinsic characteristics is the 
length L of a system appears now more pronounced. 

In order to calculate the long-range averaged characteristics of the energy spectra of 
([5]) we numerically diagonalized the Hamiltonian matrix in the representation of the base 
Fock states of the vibrons 1 and 2. Taking A^i and A^2 base Fock states for each vibron 
we introduced an appropriate ordering of sites [21j so that a state vector turned to be a 
vector with Ni ■ N2 elements. Numerical determination of high energy states produces 
an inevitable error because of truncating the base space. Varying the numbers of the 
base states we investigated the convergence of the results. For practical purposes we 
limited ourselves by the base size A^^i x A'"2 = 75 x 75 from whence about 1100-1200 states 
appeared to be trusty (giving the convergence up to O.lx (characteristic level spacing)). 
This number of energy levels was used in the following calculations. Thus, the obtained 
spectrum needs to be unfolded in order to ensure its homogeneity (that is normalized so 
that the local averaged level spacing is unity). To do so it is necessary to fit the 'staircase' 
N{E) (or the corresponding level density p{E) = J2i ^{E — Ei)) by a reasonably chosen 
smooth function. Different choices of the latter are possible. The calculations below 
used the unfolding by fitting the level density by a third-order spline polynomial (after 
smoothing the data via Gaussian smoothing with the smoothing constant 0.3) in the 
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domain of interest, that is for level numbers 100 — 1200. In this interval the level 
density of our two-vibron system grows almost linearly with the energy (we note by 
passing that this behaviour is in a qualitative accordance with the semiclassical Weyl 
formula giving generically a linear growth of the level density for a two-vibron system) , so 
that such fitting turns to be more than satisfactory. We performed sample calculations 
with an increasing order of polynomials, and using other algorithms of producing the 
smoothed level density p^m as well (in particular, varying the parameters of the Gaussian 
smoothing and using the smoothing via moving average/median); the results appeared 
indistinguishable within statistical errors. 

In order to compile a statistical ensemble out of the (deterministic) quantum 
problem one has to perform averaging of the numerical quantities of interest over 
different parts of the spectrum taken as members of this ensemble. We included into 
the statistical analysis the energy intervals centered at different starting values ranging 
between 100-1100 levels which corresponded to the energy values between (—10) and 
(—5) and between 25 and 30 depending on the values of (a,/5). However, as it follows 
from the results below, the ensemble averaged properties can show a bias (sensitivity 
to the starting point in the energy range), therefore this procedure must be followed 
with maximal care. To improve the statistics we also followed the standard procedure of 
collecting statistical data from small intervals in the space of parameters (a,/5). To be 
exact, for a given point of interest (ao, Po) we collected the energy values from the square 
(ao, Po) ± 0.25 in the parameter space with the step 0.1. (The exception of this rule was 
collecting the data for the points where ao = (3q. There, this restriction was followed 
for choosing the neighboring parameters for the sampling) . From figure 4 it follows that 
indeed the averaged statistical data present a slight bias when taking different energy 
intervals. Therefore we presented the results in figure H] separately in two energy bins for 
lower (100-600) and higher (600-1100) levels. It is worth noting that the nonsymmetric 
parameter values (a ^ (3) appear to present more robust statistics with respect to the 
energy interval than the symmetric ones (a = (3). The same conclusion was already 
drawn from the investigation of the short-range statistics [U |T7j where we showed that 
the spectral properties of the nonsymmetric JT models were more homogeneous with 
respect to the choice of the energy interval. Taking smaller energy bins than those of 
figure 4 however did not change much the appearance of the corresponding figures. 

In figure [3], the samples of the long-range spectral correlations A3 of the U- 
transformed Hamiltonian ([5]) (unfolded in a fashion described above) are presented for 
several sets of interaction parameters. It is seen that the spectral rigidity shows up 
serious deviations from the Poisson behaviour as well as from the logarithmic dependence 
expected in the frame of RMT. For large n it supposedly tends to a saturation, and a 
characteristic linear domain is evident for small n (up to n ~ 15). 

In figures m and [5l we plot the level compressibilities (slopes dK^{n) / dn\n=Q of the 
A3 curves) for a range of parameters The slopes in figures H^, |4)d were calculated 
from averaging over 500 levels in the ranges 100-600 and 600-1100, respectively. Taking 
lower levels would bring us below the 'diabatic line' where the system feels only one 
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Figure 3. Spectral rigidities (Aa-statistics) in different segments of the spectra for 
the E(X)(&i + 62) JT model. The slopes are linear for n at least up to 15 and tend to 
nonuniversal saturation values A^maxiP) for large n. 



potential well of the 'effective potential' (see [IHl [IT] for a detailed discussion). It is 
seen from figure 4 that results for different intervals of energy show differences, but, 
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Figure 4. Slopes dA3{n)/dn\n^o (level compressibilities') of the long-range statistical 
measure for different model parameters {a, (3) for different parts of the spectra. 
The grid lines indicate the Poisson line (1/15) and the marginal semi-Poisson limit 
(x = 0.5/15 ~ 0.033). The dispersion of the probability distributions pertaining to 
the displayed mean values is ~ 0.01. The pictures show a strong dependence on the 
position of the fragment within the spectra. 



nevertheless, similar tendencies. Apart from the cases close to the special symmetries 
{a = (3 and a -C /3, a S> /3) the slopes show a markable accumulation close to the value 
X = 0.5/15 (the horizontal grid line at 0.033) indicating the semi-Poisson limit. This 
limit is more pronounced at higher levels (interval 600-1100), where even the points 
corresponding to the symmetric cases a = j3 are close to this line. For lower energy 
interval (levels 100-600) this universality is seen less, although all A3 curves points 
are far from both Poisson and RMT limits. For lower energies there is also a more 
pronounced difference between symmetric and nonsymmetric models: meanwhile the 
nonsymmetric values of parameters still tend to the semi-Poisson limit, the symmetric 
points however show essential deviation from this limit towards the Poisson case. 
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In figure [3], the A3 curves for /? = 2, 4 for a = 3 (scaled hj Q = 1) averaged over 
the whole relevant extent of the spectra (up to the level n = 1100) are situated very 
close one to another, i.e. they are getting close to a limit with the minimal slope. 



dAs /dn 

. 08 
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Figure 5. The cumulative slopes dA3{n) / dn^n^Q averaged over the relevant extent of 
the spectra. Large deviation of the points a — 4, (3 > A from the suggested marginal 
line is explained by their relative proximity to the symmetric point a = /3 = 4. 



The said marginality also emerged in the statistics of NNS [T]. We have shown 
that the dispersions o"^ of NNS distributions in the range of parameters far from the 
mentioned special symmetry cases tend to a limit cr^ ^ 0.5 characteristic for the M-I 
transition of the Anderson model. 



3. Quasiclassical description in the phase space 

It is interesting to discuss the mapping between the exposed results and the semiclassical 
description of the same system. The search for the quantum chaotic patterns historically 
meant looking for quantum correspondence to the chaotic behaviour of the classical 
trajectories. First of all it is to be noted that the passage to the (semi) classical 
description in a two- or many-level electron-phonon system can be performed in several 
fashions. The most straightforward "semiclassical" description of the quantum JT 
problem is obtained if we pass to the real Bloch variables made of pseudospin Pauli 
matrices: 

x{t) := {a,)t y{t) := {ay)t z{t) := {a,)t (7) 

and perform the appropriate decoupling of the oscillatory (phonon) and electron 
variables in the initial Hamiltonian (jl]) within the framework of the Born-Oppenheimer 
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approximation: {Qid^) {Qi){t)x{t) etc. Another way of passing to the semiclassical 
description follows from the form of the Hamiltonian ([5]) where the electron degrees 
of freedom are exactly eliminated at the cost of the strongly nonlinear coupling of the 
phonon degrees of freedom. The semiclassical decoupling then means the decoupling of 
the vibrons 1 and 2. The classical equations of motions in this variant are implied by 
the following classical Hamiltonian as a function of the coordinates 7j and respective 
momenta TTf. 

Hill, TTi, 72, TTa) = I IttI + 7tI + 7i + 72I +«7iT/?72 exp (-2(7^ + vr^)) (8) 



2 

The trivial linear stability analysis shows that the stationary point (7i,72,vri = 
TT2 = 0) is unstable unless a 3> /?, a -C /3 or unless the parameters (a,/?) are located 
in the small domain near the point a = f3 = 0. The crucial characteristics of interest 
is however the chaoticity of the phase space, in particular, its fraction fi occupied by 
the chaotic trajectories. The formulae which relate this quantity to the parameters of 
the statistics of quantum levels for both short- and long-range statistics are widely 
known [SI HI [5]. We investigated the chaoticity of the trajectories implied by the 
Hamiltonian ([8]) and performed rough estimations of the 'chaoticity parameter' /i for 
different energies. We studied the sensitivity of the classical equations of motion to 
the initial conditions checking whether two initially neighboring trajectories diverge in 
the course of time (that is studying the Lyapunov index [U [5l |22]). The fraction of 
the chaotic trajectories at given energy E was then estimated by taking at random the 
initial points respecting if (71, tti, 72, 712) = E. The detailed quantitative presentation 
of this work is to be given elsewhere; for the sake of the present contribution we just 
mention the results briefly. The classical phase space in the energy range corresponding 
to the energies of quantum levels used (that is up to ~ 30) appears to be pronouncedly 
chaotic. The fraction of chaotic levels increases rapidly with the increase of energy and 
reaches its maximum /i ~ 0.90 — 0.99 for the energies typically E ^ 10 which roughly 
corresponds to the lower energy bin of figure H] (levels 100-600). With further increasing 
energy the value of fi then slowly decreases. Its estimations for E between 20-30 range 
from 0.3 to 0.75 strongly depending on the parameters a and f3. This latter energy 
interval approximately corresponds to the higher energy bin of figure |Dd. We note once 
again the more pronounced universality of the quantum characteristics at the upper bin 
(figure Hb) which appears to correspond to the moderate values of the degree of the 
classical chaoticity. The quantitative conclusions for this item is however a challenge for 
a separate study. 

In any case, the above results confirm the nonuniversal degree of the classical 
nonintegrability demonstrated by the dependence of the chaoticity parameter fi on 
the model parameters a and (3 and on the location of the energy segment in the 
spectra. Moreover, the linear dependence of the spectral rigidities in certain parts of the 
spectra imphes a strong suppression of the level repulsions because of the presence of 
regions with various degree of chaoticity and related chaos-assisted tunneling between 
the regions of small degree of chaoticity. Analogous nonuniversal behaviour of the level 
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compressibilities presented in figure H] indicates a cumulation of most of their values in 
the chaotic region (0,0.033) close to the semi-Poisson values 0.5/15 depending on the 
spectral segment where the chaoticity parameter is correspondingly reduced. 

4. Fractal dimensions and multifractality of the JT-excited wavefunctions 

An independent quantitative analysis of the properties of the £0(61+62) excited spectra 
can be obtained by exploring the scaling of the inverse participation ratio ([2]) and ([3]) 
in the spectral representation of Fock states and following calculation of related fractal 
dimensions. In figure 6, we give the samples of the scaling of log J2,L(n) (from equation 
([3])) for several successive states as functions of the box size / (L ~ 1/Z^) comprising P 
Fock states. The linear slope in the log-log coordinates stretches up to the box sizes 8x8 
which indicates the existence of a well determined quantity D2 for each level, although 
the generalized fractal dimensions show slight level-to-level fluctuations. 



a=3 [5=2 



-2 



-5 



10 


8 






5 


4 


3 


2 


1 


★ 

□ 

- □ 

♦ □ 
♦ 

♦ 


□ 


★ 

□ 


★ 


★ 


★ 








♦ 


♦ 


♦ 


□ 
♦ 


□ 
♦ 


□ 


★ 

□ 


★ 














♦ 


♦ 


□ 


















♦ 




10 


8 






5 


4 


3 


2 


1 



box size 



Figure 6. Examples of scaling of the inverse participation ratio /2^(f^) as a function 
of the box size / in the box-counting algorithm as the illustration of the reliability of 
the scaling ([3]) for our model. Different shapes of symbols correspond to the values 
for different levels n: diamonds for n — 153, stars for n = 157, squares for n = 161. 
The domain of linear slopes (proportional to the fractal dimensions D2) extends up 
to the box sizes ~ 8 x 8 (at larger box sizes the linearity is violated due to the size 
effects). Slight level-to- level fluctuations of the slopes are also seen. Similar scaling 
can be demonstrated for q = 3, 4 as well. 

In figure U\ the averaged fractal dimensions Dg, q = 1,2,3,4, in the spectral 
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representation are displayed for a set of interaction strengths. Small variations of Dg 
{Di < Dj ai i > j) testify a weak multifractality of the respective wavefunctions. The 
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Figure 7. Fractal dimensions Dg averaged over n for different pairs oi a ~ 3, /3. A 
weak multifractality is evident when comparing Dq for different q. 



fractal dimensions Dg show strong fluctuations with changing n (figure [8] so that one has 
to speak rather about their statistical distributions [23] . It is seen that this distribution 
may have a pronounced drift - there occurs a crossover between dimensions d = 1 
and d = 2 at low n. It becomes more homogeneous when n increases. In view of the 
interpretation in terms of the ID 'pseudolattice' in j-space with natural size of the order 
~ ^/n this can be understood as a tendency towards a universal distribution in the limit 
of a very large lattice size. Such a limiting behaviour of statistical characteristics of the 
distribution of fractal dimensions for large lattice sizes of the Anderson model at M-I 
transition was conjectured by Parshin et al [23] . 

Despite the level-to- level fluctuations of the fractal dimensions (figure [8]) their values 
averaged over level numbers appear to exhibit marginal universality in a similar sense 
as the universality marked in figure |H The averaged fractal dimensions D2 in figure [9] 
show a markable tendency to values 1.5 ± 0.1 apart from the case a P or a <ti /3. 



5. Conclusion remarks 



The statistical properties of the investigated JT model are closely related to its spatial 
symmetry given by the interaction strengths a, (3. Our numerical analysis of the 
spectral rigidity, its slopes and fractal dimensions bring the evidence for their apparent 
nonuniversality and spectral inhomogeneity. Namely, our results allow us to conclude 
that the nonuniversality is related by the dimensional crossover between d = 1 and d = 2 
when changing the parameters a, (3. A support for this suggestion is the behaviour of 
the fractal dimension D2 (figure [9]) which approach very close to the value 1.5 at the 
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Figure 8. Fluctuations of fractal dimensions D2 over levels n for different pairs a, /?. 

points a = f3. In the neighbourhood of this point the nonuniversahty of D2 is the most 
moderate; the deviations from 1.5 extend within the interval ±0.1. 

We have shown that although the long-range correlation measure A3(n) and the 
partial compressibilities (averaged over partial segments of the spectra) are nonuniversal 
and highly inhomogeneous, the slopes averaged over the whole relevant spectra approach 
marginally (quantitatively close) to the universal value characteristic for the semi- 
Poisson distribution, 0.5/15 ~ 0.033. This behaviour is analogous to the well defined 
marginal behaviour of the average fractal dimensions D2 of the respective wavefunctions 
mentioned above. To our knowledge, till now the theory justifying the existence of the 
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limiting 'universal' value related to the said distribution is lacking. The only relation 
X = {d — D2)/2d derived by Kravtsov and Lerner [15] was the first bridge between the 
fractal dimension D2 and the long-range level statistics. However, it is proven to be 
valid only for small values of compressibility and does not fit to our calculations. 

We have also revealed the corresponding nonuniversal and inhomogeneous 
behaviour of the chaoticity parameter - the fraction of the chaotic phase space of the 
classical trajectories. The detailed analysis of this parameter in the whole phase space 
would be very useful for further considerations of phenomena especially related to the 
mixing [7] of the phase space between the regions of different degrees of chaoticity. The 
presence of the mixing phenomena is supported by the multifractal behaviour of the 
wavefunctions. Because of the nonexistence of the extended states in our model we 
can conclude on the concept of the chaos-assisted tunnelling between the domains with 
reduced stochasticity through those with high degree of stochasticity. 

Besides the Anderson model close to the metal-insulator transition [8j another 
models sharing the universal semi-Poisson statistics are known - in particular, the 
Bogomolny [10] plasma model with screened Coulomb interactions and the 'daisy 
models' [20]. 
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